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Introduction 


This  paper  deals  with  a  cethod  of  calculating  th*.  «„i  the 

path  of  an  orbiting  body  fron  a  nocinal  or  reference  trajectory.  The 
fero  in  which  the  solution  is  cast  was  motivated  by  .*  f-rsicalar 
perturbation  problem.  Stanford  University  is  developing  a  '"drag-free", 
or  ‘‘drag-makeup”,  scientific  satellite  which  is  designed  to  follow  a 
purely  gravitational  orbit.'*'  The  satellite  consists  actually  of  two 
satellites:  an  inner  sphere  or  proof  mass,  and  an  outer  concentric 
shell.  The  relative  position  of  the  shell  with  respect  to  the  inner 
sphere  is  sensed  with  a  capacitive  pickoff .  The  position  signals 
command  an  active  translation  control  systea  which  fires  jets  Eoanted  on 
the  outer  shell  so  that  it  chases  the  inner  sphere  without  ever  touching 
it.  Thus  the  proof  nass  is  shielded  froc  gas  drag  and.  solar  radiation 
pressure  and,  except  for  very  sssall  disturbances  caused  by  force 
interactions  with  the  outer  shell,  it  follows  a  purely  gravitational 
orbit. 

The  problem  which  motivated  the  present  study  was  tc  determine  the 

_10  -S  . 

effect  of  these  snail  disturbances  (about  10  to  10  g  }  over  tine 
periods  up  to  a  year.  Furthermore,  the  answer  was  desired  directly  in 
terns  of  the  deviation  of  the  satellite's  path  froc  the  path  which  would 
be  followed  by  an  earth  satellite  acted  upon  by  gravity  only.  Therefore, 
the  technique  of  perturbation  of  the  coordinates  was  selected  as  the  basis 
of  our  approach. 

The  technique  of  coordinate  perturbation,  which  began  with  the 
2  3 

work  of  Encke  and  Hill  in  the  last  century,  has  found  increasing 
use  in  modern  tines  for  orbital  theory.  The  linearized  perturbation 
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equati'f^s  about  a  circular  orbit  (which  arc  merely  li-.il’ s  lunar  equations' 

without  the  mutual  gravitational  teres  (see  equations  (1),  (2),  (3)  with 

*5  5*6 

e  =  0}  have  been  applied  in  recent  years  by  Wheeloa,  Geylicg,  and 

7 

Clohessy  and  Wiltshire  to  a  sssber  of  satellite  perturbation  problems. 

8  _  S.IO.U 

Battle  and  Darby  *  give  state  transition  matrices  fer  general  conic 

sections  »hich  also  say  be  applied  to  satellite  perturbation  and  guidance 

12 

problems,  and  recently  Tscfcaurer  and  Hesgpel  have  applied  the  linearized 
Hill’s  lunar  equations  to  the  minimum  fuel  rendezvous  problem. 

Is  some  types  of  orbital  problems  (as  in  the  mentioned  example  of 
determining  the  effect  of  internal  force  errors  on  the  orbit  of  a  drag- 
free  satellite),  it  is  desirable  to  compute  the  perturbations  of  the 
coordinates  when  the  satellite  is  subjected  to  very  small  disturbances 
for  many  thousands  of  revolutions.  In  this  case,  the  linearized  Hill’s 
equations  are  useful  only  for  very  very  small  eccentricities;  variation 
of  parameter  techniques  do  cot  yield  an  answer  directly  in  the  desired 
fora  (i.e. ,  as  deviations  of  the  coordinates);  and  direct  numerical 
integration  proves  both  costly  and  inaccurate,  when  carried  out  over 
long  time  intervals.  Hence  a  different  approach  is  sought. 
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The  Tschauner  and  Herrpel  Equations 

13 

Tschaucer  aed  Kecpel  have  shovs  chat  if  the  &<r jliz*.  f  r  it 
equations  of  cotx^n  are  linearized  about  a  ry^ital  elliptical  ■  r_.il  m 
a  rotating  reference  frace  (see  Appendix  A) ,  the\  assu-ic*  the  .vry 
simple  fora: 


m*9 

5 


i  e  cos  0 


§  -  2  /  =  cr 


(1) 


where  § 


P  ,  P  ,  P  are  snail  perturbing:  accelerations  along  the  u  > 
12  3  I 

u  ,  a  axes  respectively, 

*  3 

B  is  the  instantaneoas  radius  of  the  nominal  elliptical  orbit  , 

0  is  the  true  anomaly  in  the  nominal  orbit, 
e  is  the  eccentricity  of  the  nominal  orbit, 

=  ©»  the  tine  rate  of  change  of  true  anomaly, 
u  »  a  c  are  relative  coordiantes  shewn  in  Figure  1,  and 

**  £  3 

the  prime  (')  signifies  |  - 

2  2  2 

In  deriving  these  equations,  terms  of  order  -  ,  -  .  ^  and  higher  are 
neglected.  If  the  equations  of  motion  in  cylindrical  form  are  linearized 
as  shown  |s  Figure  2,  with  5  =  |,  %  =  c„  and  £  as  before,  equations 
(1)  through  (3)  are  again  obtained.  Sow,  however,  -  say  be  arbitrarily 
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r 


2  2  2 

large  vhile  t eras  of  order  «  ,  -J  .  £  ,  and  higher  are  neglected. 
Equation  (3) ,  of  course ,  represents  sisple  out-of-plane  hanaonic 
notion  and  needs  no  discussion. 
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Solution  of  the  Yschauner‘>3empel  Equations 


By  introducing  matrix  notation  and  defining  the  system  state  matrix 
x(B)  to  be 


x(a)  = 


(4) 


equations  (1)  and  (2)  may  be  combined  and  written 


x*  (9)  =  F(0)  x(9)  +  B (9)  u (9) 


(5) 


where 


<6> 


(V 
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■ 


•*: 


and 


»>  ■  r-”) 

vp2W 


(8) 
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It  is  well  known  from  the  theory  of  Floquet  (see  Appendix  B)  that 

a  system  governed  by  equation  (5)  where  F(9)  =  F($  +  2:t),  has  a  state 

* 

transition  matrix,  X(g ,  9Q),  which  can  be  written  as: 

x(9,  eo)  =  R(e,  0o)  eB(e  "  eo)  (9) 

where  R (9,  t?Q)  =  R (6  +  2it,  9^)  is  a  periodic  4  *  4  itrix,  and 

B  =  in  X(g  -t-  2k,  9  )  is  a  constant  4x4  matrix 
2jt  o  o 

whose  eigenvalues  determine  the  system  stability. 

The  unforced  part  of  equation  (5)  is  said  to  be  kinematically  similar 
to  the  constant  system 


w’  =  Bw  . 


(10) 


F<9),  R(9,  e  )  and  B  are  related  by 
o 


B  =  R  1(9,  6o)F(9)R(9,  0q)  -  R  1(0,  90)R'(9,  6q)  , 


(11) 


and  equations  (10)  and  (11)  are  known  as  the  Lyapunov  reduction  of 
equation  (5).  By  an  appropriate  linear  constant  transformation 

z  =  Qw  (12) 

♦Formally,  the  state  transition  matrix  of  an  n^-order  linear  system  of 
differential  equations  in  first-order  matrix  form  is  an  n  x  n  matrix  whose 
columns  are  n  linearly  independent  solutions  of  the  free  equation,  such 
that  X’(9,  90)  =  F(9)X(9,  9Q)  and  X(9Q,  $Q)  =  U,  the  unit  or  identity 
matrix  (see  Appendix  B). 
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(where  Q  is  a  constant  4x4  matrix) 

equation  (10)  may  be  transformed  into  its  Jordan  normal  form: 

z*  =  A 2  (13) 

where 

A  =  QBQ_1  (14) 

The  eigenvalues  of  A,  together  with  the  structure  of  the  Jordan  blocks 
determine  the  stability  of  the  free  solution  (u(0)  =  0)  of  equation  (5), 
and  it  is  possible  to  give  the  state  transition  matrix,  X(9,  0Q)» 
directly  in  terms  of  A: 

Xfe,  e  )  =  P (0)eAfe  "  9o>P~1(9  )  (15) 

o  o 

where  x(9)  =  P(9)z(9).  (See  Appendix  B)„ 

The  periodic  part  of  the  state  transition  matrix,  R (9,  0q)  is  given  by 

R (9,  0  )  =  P(9)P-1(0  ),  (16) 

o  o 

and  furthermore 

Q  =  P_1  (9  )  .  (17) 

o 

13 

It  has  been  shown  by  Tschauner  and  Heap el  (who  have  obtained  the 

matrix  P  *(9)  in  closed  fora),  and  also  by  the  present  authors,  that 
equation  (10)  is  kinematically  similar  to  equations  (1)  and  (2)  with  the 
Jordan  canonical  form  of  B  given  by 
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♦The  noroal  fora  (equation  (13)  with  A  given  by  equation  (18))  corresponds 
to  two  decoupled  second-order  systens:  a  pure  inertia  or  1/s2  plant  and 
£  I  araonic  oscillator  with  a  natural  period  equal  to  the  orbit  period.  The 
1^  2  plant  »ay  be  interpreted  physically  as  motion  in  a  similar  coplanar, 
coaxial  ellipse  with  higher  or  lower  total  energy.  The  harmonic  oscillator 
corresponds  to  motion  in  a  coplanar  ellipse  with  the  sane  period,  but 
with  different  eccentricity  and/or  orientation. 

♦♦This  remarkable  property  is  not  usually  possessed  by  even  the  simplest 

of  periodic  systems.  Compare,  for  example,  Mathieu's  Equation, 

0  +  o  2(1  -  e  cos  2  ut)  0=0. 
o 


-  9  - 


The  aatrix  P  (9)  is  given  by  Tschauner  and  Heapel: 


-2q,  4  ep‘ 


e  sin  0  -£(1  +  e  cos  0)  £e  cos  0 


=|(3  4-  e  cos  0) 


e  sin  0  -$(2  4  e  cos  0)1 


where: 


c  =  ^[.1  -  (1  4  2e2)i|l-e2  ’sin  Q  -  (2  4  3e  cos  @  4  e2)sin  (23) 


_ 

P2  =“ ~U  4  sjl^e2)  -  -  Cl  -  e2)  ]cos  0 


4  \  [(1  4  2e2)fl~e2  -  ljcos  2  0  -  e  p  sin  1 
o 

=  1(2  4  e2)^l-e2  -  2]sin  Q  4  jr  [(1  4  2e2)^|l-«?  -I’sin20 


,,  „\2  -  -1 

4  (1  4  e  cos  0)  sin 


qj  =  (1  4  a  cos  ©)  , 

pi  =  sin  0(1  4  e  cos  ©), 

_  sin  0  [e  4  (1  -  COS  0  ] 

^  -  1  4  e  cos  0 


If  0  is  chosen  to  be  zero,  then 
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/ 


p'1  <0)  = 


-(1  e)  (2  +  e) 


-J(3  *  e) 


(e-f-1)  r-l+(l-e)2\l-e2l 


3e 


-1(1  -5-  e) 


\ 


3 


0  -(1  +  )' 


¥ 


0  -|(2  4-  e) 


(29) 


/ 


/ 


P(0)  = 


3e 


2  3/2 

(1-e)  (1-e  > 


(r-e)2^ll-e2 


2  ■»  e 
1  +  e 


3  -fr  e 
1  +  e 


0 


-2 


(1-e)  (l-e2)3/2 


-2>  l-(l-e)24l-e2] 
e (1-e) 2^l-e2 


\ 


2(1  -r  e) 


-2(2  4-  e) 


(30) 


Fro®  equations  (14)  and  (17) 


B(e)  =  P(0)A  P_1  O'3) 


(31) 


so  that 


-  11  - 


and  also  B  =  B(0)» 
o 

It  aay  be  seen  by  direct  differentiation  that  the  solution  to 
equation  (5)  is 


x(9)  =  X(9,  9  >x(9  )  +  Xfe,  e  ) 
o  o  o 


X-1(t,  9  )D(t)u(t )dT 
6 


C33> 


where  is  the  initial  value  of  the  system  state  matrix  and 

X<9,  eo)  =  Rfe,  0o)€B(9“9o)=  PfcJe'^oV1^)  .  C34) 

(This  solution  aay  also  be  obtained  by  variation  of  paraceters.  See 
Appendix  B)  If  one  attempts  to  use  equation  (32)  directly  to  determine 
the  effect  of  saall  perturbing  accelerations  over  nany  revolutions, 
serious  numerical  difficulties  are  encountered  which  result  both  in  loss 
of  accuracy  and  in  excessive  computation  tine. 
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Solution  for  Constant,  Periodic,  and  Algos t-Periodic  Perturbing 
Accelerations 

As  in  the  case  of  computing  perturbations  for  a  drag-free  satellite, 
it  often  happens  that  perturbing  accelerations  are  constant  or  periodic. 
It  is  then  possible  to  compute  their  effect  at  any  future  tine  merely 
by  computing  their  effect  over  one  orbit  revolution.  If  the  disturbing 
acceleration  has  the  fora 


u£©)  =  uC 9  -5-  2s) 


(35) 


it  can  be  shown  (see  Appendix  C)  that  the  solution  to  equation  (5/ 
[that  is,  equation  (33)]  can  be  written 


X 


xfe)  =  X(©  -2ar  X, 


©  )£*  x<©  )  *  x<©  -  2a  x,  ©  ) 

o  o  o 


(£') 

k=l 


-s-  x<a  -  2x  x,  ©  ) 

O 


©  )B(t)a(t)dt 
o 


(36) 


where  X  is  the  largest  cumber  of  complete  revolutions  in  (9  —  ©  ), 

o 


C  =  Xfe  s-  2s,  ©  ),  and 

o  o 


©  -52* 


a 


(t,  ©  )E(T)u(r)d-- 
o 


(37) 


The  solution  as  gives  by  equations  (36)  acd  (37)  requires  integration 

over  a  maximum  of  one  orbit  revolution,  regardless  of  the  actual  number 

of  orbit  revolutions  contained  ic  the  range  of  interest  (©-©)-  Thus 

o 

the  difficulties  mentioned  in  the  application  of  the  solution  in  the  fern 
given  by  equation  (33)  are  overcome.  The  restriction  of  the  disturbance 
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to  constant  or  periodic  in  Q  case  can  be  relaxed 


Ii 


U (9)  =  ufe  -5-  is  JO 


(38) 


where  X  is  an  integer,  it  can  be  shown  (see  Appendix  €)  that  the 
solution  [equation  (33)  ]  can  be  written 


xfe)  =  X<3  -2s  X,  0  )C  x(5  )  f  X(9-  2s  X, 

o  o 


r-l 

? 

k=0 


*  X(S  -2s  X, 


e  )c 

C 


S-ri£ 


/ 


£-2srX 


X-1(S, 


g  )B(?)a(t)dt 
o 


(39) 


where  X  is  the  largest  crusher  of  complete  revolutions  in  (g  -  ©  ), 

o 

r  is  the  largest  integer  <  5/K, 

C  =  X<9  -s-  2s,  g  ),  and 

o  o 

/o*2sX 

X_1(s,  ©o)B(x)a(s)dr  (43) 


In  this  case  the  solution  over  anv  interval  (9  —  0  )  requires  integration 

o 

over  a  maximum  of  i  revolutions.  Thus  the  constant  r  defined  above 
is  a  figure-of -merit  for  the  solution  in  this  form.  The  larger  r,  the 
more  relative  value  equation  (39)  has  over  equation  (33). 

(he  farther  generalization  of  the  forms  of  the  perturbing  acceleration 
can  be  made.  If  instead  of  equations  (33)  or  (38)  we  have 


afe)  =  u(9  -5-  6)  (41) 

where  8  f  2x  K  for  X  =  0,  1,  2,  ....,  the  solution  may  be  approximated 
(again  see  Appendix  C)  as  closely  as  desired  by  selecting  an  integer  S 
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such  that 


K  t  =  2fi  M  (42) 

for  socse  integer  M.  Then  the  solution  to  equation  (5)  is  again 
equation  (39),  with  N,  C,  and  I  as  defined,  but  with  r  an  integer 
such  that 

r  K  fc  <  2s  X  <  (r-rl)  K  &  (43) 

In  this  final  case  integration  is  required  over  a  maximum  of  M 
revolutions*  Of  course,  the  larger  the  selected  value  of  S,  the 
greater  the  accuracy  obtained  in  the  approximation  of  equation  (42). 

The  usefulness  of  the  solution,  in  this  case,  is  dependent  upon  the 
nature  of  the  actual  problem. 
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Restriction  of  Initial  True  Anocaly  to  Zero 


If  the  initial  value  of  the  true  anocaly  is  taken  »•>  r  i  .  -  *i| . 

no  real  restriction  of  the  general  problem  is  Thi  -  i-  s,  L-t-cause 

stipulation  of  9q  =  0  sicply  requires  a  compensatory  adjustment  in  the 

initial  value  of  the  system  state  matrix  x(v  ).  Then  the  solution  for 

o 

perturbing  accelerations  of  the  form 

u (.6)  =  u( 6  4-  2-r)  (35) 


can  be  written  in  a  manner  especially  adapted  for  rapid,  accurate  evalu¬ 


ation.  If  6=0  and  x(6  )  =  x  ,  equations  (36)  and  (37)  become 
o  o  o 


x(6)  =  X(j,0)Xk  (2^,0)x  +  X(a 

o 


/v  *  \ 

,0)J  2,  X  <2^.°>} 

\k=l  / 


f°  -1 

X(o,0)  /  X  (t,0)D(T)u(:)d 


(44) 


where  N  is  the  largest  number  of  coatplete  revolutions  in  0, 


a=0-  2j3f 


(45) 


.2-r 


rl, 


1=1  X  (T,0)D(t)u(t)dt 
y0 


■X 

/  i 

-fce  (2+e) 
(1-e)  2Xl-e2 


X(2j,0)  = 


(l-e)2Vl-e2 


1  0 


0  0 


(46) 


\ 


-6-ne  (I-~e) 


(1-e) 


-6?  (2-te)  (1-te)  Q  -6j(l*e)2 


(l-e)2^l-e2 


(47) 
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X~* (c,0)  = 


(1,1) 

(1,2) 

(1,3) 

(1,4) 

(2,1) 

(2,2) 

(2,3) 

(2,4) 

(3,1) 

(3,2) 

(3,3) 

(3,4) 

(4,1) 

(4,2) 

(4,3) 

(4,4) 

where 


(1,1)  = 


4  -f  e  -  3cos  g 
1  +  e 


,,  os  “sin  o  (1  +  e  cos  a)  , 
11,41  ”  1  e 

(1,3)  =  (2,3)  =  (4,3)  =  0  , 

,,  2  +  e  -  cos  a(2  +  e  cos  a) 

a'4)  *  - r~Te - 


(2,1)  = 


(2,2) 


(2,4)  = 


3e(2+3e  cos  g  +  e2)  (cr-sin  *X)  -  (3+6e2)4l-e2 


sin  o  , 


u,  2.3/2 
(1  -  e)  (1  -  e  ) 


2  _  j  r,,-“  ■  3  2  3  2 

3e _ sin  rrd+e  cos  a)  (o-sin  X)-?!-^  [  2e+e  -(1-e  )cos  q-(e-«-2e  )cos 


a  -  exi  -  =2,3/2 


cos  g| 
(54) 


2  —1  ( - 2,  2  3 

3e(l-K?  cos  g)  (q-sin  \)-11-e  [  (2+e  )sin  p+(e*2e  )sin  g  cos  a]  , 


(1  -  e)  (1  -  e2)3/2 


(  .  _  3(2+3e  cos  g+e2)  (q-sin  3X)-(6+3e)Tl-e2  sin  g 

-  2! - 0  » 

( i  .  Ji  _ 


(3,2)  = 


(1  ~  e)  •fl  -  e^ 

_1  j  2  2  2  2 

3e  sin  o  (1+e  cos  j)  ( i-sin  \)-ll-e  [2+e  -(2-2e)cos  q-(2e-fe  )cos  g]  , 

(1  -  e)2i|l  -  e2 


(3,3)  =  1  , 
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(3,4)  = 


(4,1)  = 


3{l+e  sos  cf)^(o-sin  (4-e)sin  rt(2e'H?‘')sin  ~  c< 


(1  -  e)  1  - 


S(cos  t  -  I) 
1  +  e 


(59) 

(60) 


(4,2)  = 


2sin  g(l  +  e  cos  i? 
1  +  e 


(61) 


(4,4)  =■ 


-(3  +  e)  -s-  2cos  g(2  e  cos  >7) 
1  +  e 


(62) 


ana 


_  sin  gfe  -?•  (1  -  e~*)cos  q] 

*v  ~  1  +  e  COS  <y 


(63) 


If  J  is  the  Jordan  canonical  fora  of  X(2a,0)  (given  by  equation  (47)  ) , 
then  J  is  jiven  by 


1  . 


J  s 


0 


2jt 


0  0 


0 


0 


(64) 


and  (see  Appendix  3) 


X(2jt,  0)  =  P(0)JP_1(0) 


(65) 


-1 


where  P(0)  and  P  (0)  are  given  by  equations  (29)  and  (30) .  Noting 
then  that 


(2 


a,0)  =  P(0)jV1(0) 


(66) 
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and  defining 


-l* 


(67) 


equation  (44)  can  be  written  in  a  form  which  is  convenient  for  calculating 
x($)  when  N  is  large: 


x(e)  =  x(a,o)p(o)jNp  x(o)  +  x(o,0)p(0)sp“1(0)i 


+  X(a,0) 


(t. 


0)D(r)u(T)dT 


(70) 


where  N  is  the  largest  number  of  complete  revolutions  in  fi, 
X  1(t,0)  is  obtained  from  equations  (48)  through  (63), 
D(t)  is  given  by  equation  (7)  , 
u(t)  is  given  by  equation  (35), 
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(71) 


(72) 


Sampled-Data  Solution 

If  in  the  general  solution  (equation  (70))  is  re.-trictcd  to 
zero,  an  expression  is  obtained  which  represents  sanpled  value-  of 
the  perturbed  motion  taken  at  intervals  of  2x: 

x(2jt  X)  =  P(0)JN  P_1(0)x  +  P(0)  S  P_1(0)I 

o 


where,  from  equations  (29),  (30),  (71)  and  (72)  we  obtain 


P(0)J*P  x(o) 


0  0 


6en(2+e)N 

(l-e)24i-e^ 


6rr(l-fre)  (2+e)X 
(l-e)2| 1-e2 


0  0 


6eg(l*e)S 

(l-e>2{l-e^ 


S^(l-fe)  y 

(l-e)2jt? 


0  0 


and 


/ 


P(0)SP  2(0)  = 


N 


0  0 


X 


3ex(2+e)(X+l)N  „  „  3ex(l+eHX+2)X 

£%  U  — 


(1-e)2^ 


3x(l+e)  (2+e)(X+l)N  3K(l+e}2(X*-l)N 

v  ii 


(l-e)2{1^2 


(l-e)2^l-e^ 


0  0  0 


X 


/ 


(73) 


(74) 


,  (75) 
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and  where  I  is  defined  by  equation  (46), 

xq  -  x(0)  the  initial  value  of  the  systea  state  natrix. 

The  I-aatrix  has  been  evaluated  (prinarily  by  contour  integration) 
for  the  case  of  accelerations  constant  in  the  rotating  reference  frane; 
that  is,  for  accelerations  of  the  fora 


’a  » 


u(6)  = 


(a. ,  a  constant) 
5  rl 


(76) 


The  result  of  this  evaluation  is  contained  in  Appendix  D.  Using  this 
result  one  obtains  for  accelerations  described  by  equation  (76): 


/ 


P(0)SP~3  (0)1  =  — 

k/p 


(4-e)?-X 


\ 


(l-e)(l-e2)3/2  1 


3e(e~2e^2)nX 


2  2 
6es.\ 


,,  2.5/2  £  ,,  2-2  n 

(l-e)(l-e  )  s  (1-e) (1-e  )  1 


(e  ~10ef4)^X 

r,  -2,.  2-3/2  a| 

(1-e)  (1-e  )  51 


(1-e)  (1-e  ) 


2X  2r, 


3(e  -2e-2)gX 

a-o2)5/2  aT 


(7  7) 


/ 


where  k  is  the  gravitational  field  constant 

p  is  the  seailatlsrectua  of  the  noainal  ellipse 


A  closed-fom  of  the  I-natrix  has  also  been  obtained  for 
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accelerations  of  the  fora 


u(0)  = 


„  J'V 

v-  1 


V  2 


(78) 


where  K  and  K  are  integers, 

A  xS 

c.  and  c  are  arbitrary  coaplex  constants, 

A  Z 

This  result  is  also  presented  in  Appendix  D.  It  is  possible  then, 
using  the  two  foras  of  the  I-natrix  (for  constant  and  for  periodic 
accelerations),  to  derive  an  appropriate  I-natrix  for  any  disturbance 
which  can  be  expanded  in  a  Fourier  series  in  true  anocaly. 
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FIG.  3.  ACCELERATION  CONSTANT  IN  ROTATING  REFERENCE 


Application  to  the  Drag-Free  Satellite 


As  cent ioned  in  the  Int  redaction,  the  proof -rass  •  r  inn-  r  -7;-*  r 
portion  of  the  drag-free  satellite  follows  a  pare  grav.it>  ».rhi:  except 
for  very  s~all  perturbations  caused  by  force  interactions  bet seen  the 
inner  and  outer  satellites.  The  cajoritv  of  these  force  interactions 
are  essentially  fixed  within  the  satellite.  Shen  the  satellite  is  cain- 
tained  in  a  locally— level  orientation  these  forces  are  then  fixed  in  the 
rotating  reference  f raise  and  can  be  described  by  equation  (76),  i.e.. 


ate) 


(a  ,  a  constant) 
5  1 


(76) 


Plots  of  typical  perturbed  notions,  assusing  zero  initial  conditions,  re¬ 
sulting  fron  accelerations  of  this  type  (for  selected  values  of  nominal 
orbit  eccentricity)  are  presented  in  Figures  13)  through  (6) .  It  is  in¬ 
teresting  to  compare  these  plots  with  Figures  (4-5)  through  (4-8)  of 
Reference  15,  which  represent  the  solutions  for  zero  eccentricity.  As 
would  be  expected,  the  results  for  e  =  .01  are  almost  identical  to  those 
for  e  =  0,  bat  do  exhibit  the  trend  or  distortion  shown  amplified  in 
the  plots  for  e  =  .1.  The  effect  of  the  secular  terns  of  the  solution 
are  cost  easily  obtained  through  the  saspled-data  solution  of  equation  (73)  . 
Again  ignoring  initial  condition  effects,  and  selecting  for  example 
e  =  .01,  then,  using  equation  (77),  one  has 


|(2r.S) 


4rfJ 


(79) 
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STANY  IN  ROYAYXNO  REFERENCE 


IN  ROTATING  REFERENCE  F 


2  2 

-  4r.\  6-  N 

n<2-N)  - - 2% - 2*- 

k/p  ^  k  p  1 


—10  2 

If  a  =  a  =  10  ie/ sec  ,  p  =  100  dies,  plus  the  radius  of  the  Earth, 

5  n 

X  =  SOCOrev  F  1  rear) ,  and  using  the  basic  relationships  x  =  §R,  y  =  tjR, 
it  is  seen  that 


x(l  year)  =  6  n  =  20  feet 


v Cl  rear)  =  —10  ra  =  -60  dies 


These  results  verify  those  of  page  (126)  of  Reference  15. 

The  drag— free  satellite  cay  also  be  oriented  so  that  it  "-aintains 
its  orientation  with  respect  to  inertial  space.  Then  the  perturbing  ac¬ 
celeration  would  be  essentially  fixed  in  inertial  space.  If  it  is  resolved 
into  a  component  a  Iving  along  the  line  of  apsides  of  the  nominal 

orbit  and  positive  outward  (away  from  the  focus)  and  a  component  a  ^ 

* 

perpendicular  to  a  in  the  plane  of  the  nominal  orbit,  and  positive 

5 ° 

in  the  direction  of  notion,  then  the  acceleration  vector  becomes 


uCs)  = 


a  cos  Q  ~  a  sin  © 


sin  Q  +  a  cos  Q 

P 


where  a  and  a  are  constant.  Examples  of  typical  notion  are  pre- 

i°  IP 

seated  in  Figures  (7)  through  (10) .  Again  it  is  interesting  to  compare 
these  results  with  those  for  zero  eccentricity  contained  in  Reference  15 


As  was  noted  previously,  out  of  plane  notion  is  sicple,  decoupled. 


harmonic  motion,  requiring  no  discussion. 
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ACCELERATION  CONSTANT  IN  ROTATING  REFERENCE  FRAME 


as  Figures  (4-9)  and  (4-10) . 


-  30 


0.01 


where  the  acceleration  vector  a  has  been  resol .ed  as  was  done  previously 
for  accelerations  fixed  in  inertial  space . 
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FIG.  10.  ACCELERATION  CONSTANT  IN  INERTIAL 


To  determine 


u(«)  = 


H  ■  • 

\P2<3>  J 


(86) 


P1(9)  = 


It  is  possible  to  expand  u(g)  as  a  Fourier  series.  If  this  is  done  then 

1 


a,  (sin  Q .  -  sin  e  )  -  a  (cos  0.  -  cos  @  ); 


50  i  o  ^o  "1  o  | 


sin  2 g_  -  sin  29  cos  2g.  -  cos  2g 

+  0.  -  9  )  +  - 2^- - -  -  a  - =-= - 

i  o  jO  2  'o  2 


1 

i-L  (2, 

2s  1  jo 

x—  \a  (2rr  4-  0.  - 

2*  ^  rko  1 

1  V1  (  Ain  (n-s-l)  8^  ~  si 

2s  Z  |  a£°l  (n  +  1) 

n=2  1  \ 


sin  20.  -  sin  2&  cos  2 o  -  cos  2a 

.*  10  i 

%>  - »  - - - =* 


^0 


"I 

- 


.*1  )S.  -  sin  (n-f-i)0o  sin  (n-Dg^  -  sin  (n-l)©^ 

f  (n  -  1) 


(cos  (n+l)0.  -  cos  (n+l)0  cos  (n-l)0_  -  cos  (n-l)g 
j  2.  o  x  o',  . 

a  | - ; - - - ; - — - 1  >  cos  n0 


(n  -5-  1) 


x 

.  -A  {  /cos  (n+l)0.  -  CO 

-  k  1  V - 5Ti) 

n=2  1  \ 


(n  -  1) 


-  cos  (nfl)0  cos  (n-l)0.  -  cos(n-l)g 
_ 2.  * _  1  °1 


(n  -  1) 


+  a 


rjo 


/sin  (n+l)0.  -  sin  (n+l)g 

I  1  0 

sin  (n-l)0.  -  sin  (n-l)0^| 

1  (n  -  1) 

ps 

pH 

1 

a 

SjH 

1 

sin  n 9 
(87) 


COS  0 


sin  0 
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P2<9) 


(sin  0.  -  sin  Q  )  +  (cos  8.  -  cos  Q 


l-L 

2s  |  7 

^7  Sa_o(2'  + 

2s  |  *;o  i 

k  - 


«’l 


sm  20.  -  sxn  2%  cos  2a.  -  cos  20 

.  .  1  G  i  O 

9  )  +  a  - 5 - -*■  a.  - = - 

o  wfl  2  |o  2 


sin  2g.  -  sin  2s  cos  29  -  cos  2a 

V2*  +  ®i  "  *  V - ~ - -  "  “'.o - 2 - ° 


cos  @ 


sin  a 


x  ~  r  ^in  (n^-l)ai  -  si 

2s  c-i  I  ^ol  (n  +  1) 

n=2  1  \ 


I.  -  sin  (n4l)g  sin  (n- 
o 

-  — 


50 


fcos  (n+l)p  -  cos  (n-t-l)g  cos 
i  o 

(n  +  1) 


j  ^  f  /cos  (n*l)ei  -  cos  (n-r 

2s  2->  j^ol  (n  4-  1) 

n=2  1  V 


(n-l)a.  -  sin  (n-I)a\ 

— chro - j 

(n-l)a.  -  cos  (n-l)a  \) 

- - -Jj  cos  “9 

(n-l)@t  -  cos  (n-l)eA 

(a  -  1)  I 


-rl)s  cos  (n- 
o 


—  a 


(sin  (nfl)0.  -  sin  (nfl)0  sin  (n- 

2.  O 


(n  4-  1) 


(n-l)a^  -  sin  (n-l)e  \  ) 
(n  -  1)  jj 


sin  n0 

(88) 


Figure  (12)  is  a  plot  of  perturbed  notion  over  4  orbit  periods  under 
the  following  conditions: 

a  =  -a  (the  sun  lies  along  the  line  of  apsides) 

1° 

e  =  0.01 

ei  =  1350 

0o  =  2250 

In  the  nuaeric  integration  the  Fourier  expansion  was  carried  out  to  the 
19th  tern  (n  =  19) .  It  should  be  noted  that  the  paraceters  selected 
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were  chosen  Merely  to  provide  an  idea  of  the  nature  of  the  solution, 
rather  than  to  describe  sene  actual  orbit  condition.  The  probloa  of 
calculating  actual  shadow-entry  and  exit  angles  is  discussed  in  the  lit¬ 
erature  Ccf.  heference  16)  and  is  not  within  the  scope  of  this  paper. 
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FIG.  12.  EXAMPLE  OF  SOLAR  RADIATION  PERTURBATION  EFFECT. 


Application  to  Inertial  (kiidacce 


The  basic  relationship  of  inertial  guidance  is  that  geosetric 
acceleration  is  equal  to  the  output  froa  an  ideal  acceleroeeter  plus 
gravitational  mass  attraction.  That  is 


r  =  t  *  g  (89) 


vfaere  r  is  the  position  vector  of  the  vehicle, 

£  is  the  output  of  an  ideal  acceleroeeter  cn  board  the  vehicle, 

— i 

g  is  the  gravitationar  sass  attraction  vector,  and 
overscript  (r)  signifies  d/dt  in  an  inertial  frase. 

An  inertial  guidance  svstecs  coc outer  is  mechanized  such  that  it  obtains 
the  solution  to  equation  (89)  by  solving  the  Ideal  Mechanization  Equations, 
(90)  and  (91). 


c 

-» 

V 


g  -  tsx 


(90) 


c 
-» 
r  = 


u  X  r 


(91) 


-»u  “» 

There  v  =  r  , 

cs  is  the  angular  velocity  of  the  computer  frase  with  respect  to 
inertial  space,  and 

overscript  (C)  signifies  d/dt  in  the  coaputer  frase. 

20 

It  has  been  shown  elsewhere  that  free  these  three  basic  equations,  by 
perturbation  analysis,  one  obtains  the  Platform  Misalignment  Error 
Equation  (92)  and  the  Position  and  Velocity  Error  Equation  (93)  for  an 
Inertial  Navigation  Systea  in  elliptical  orbit 
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-4  -4 

w  -  e 


(92) 


=  -  K 


g 


9  ^  < 

Sr  +  oTSr  -  3u> 

5  S 


2  r 


-4 

Sr 


r  = 


“4  “■)  *- > 

-  ^  X  f  +  K 


g 


-4 
— > 
+  K 

v 


— >  — >  —4  -4  —4 

f  +  b  *{-  T)  —  25W  X  v 


d1  ^ 

v  + « (KP  ’ r  j '  x 


(93) 


where  y 


t 

g 

-4 

e 


Sr 

W  =(k/r3)^ 
s 

k 


b 


n 

— $ 
6  w 


is  the  vector  approximating  the  small  angle  rich  rotates 
computer  into  platform  axes, 

is  the  diad  representing  stabilization  gyro  scale  factor  error, 

is  the  stabilization  gyro  drift  rate  or  bias  error  vector, 

is  the  first  order  approximation  to  the  error  in  r, 

is  the  Shuler  frequency  corresponding  to  |jrj , 

is  the  universal  gravitational  field  constant, 

represents  accelerometer  bias, 

represents  random  accelerometer  errors, 

is  the  difference  between  computer  angular  rate  and  platform 
angular  i*ate  (to  first  order), 

is  the  diad  representing  first  integrator  scale  factor  error, 
is  the  diad  representing  second  integrator  scale  factor  error, 
and 

signifies  differentiation  in  an  inertial  frame. 


Comparison  of  equation  (93)  with  equation  (A10)  of  Appendix  A  reveals 
the  interesting  fact  that  the  homogeneous  form  of  the  Position  and 
Velocity  Error  Equation  of  an  Inertial  Navigation  System  in  ellipti *al 
orbit  is  identical  to  the  homogeneous  form  of  t  Basic  Perturbation 
Equation  linearized  about  an  elliptical  orbit.  It  follows  then  that 
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mc&maetn  •  jl-  ~*=v  - 


I 


equation  (93)  can  be  transformed  to  the  Tschauner-Hempel  Equations. 
If 


(94) 


coordinatized  in  a  locally-level  reference  frame,  and 

a  a  r 

6x  =  r  5|;  8y  =  r  5*1;  8z  =  r  5£;  (95) 

then  equation  (93)  becomes  equations  (96)  through  (98): 


6£"  -  1  +  e3cos  8 

-26  t)'  =  a 

(96) 

25|  *  +  6Tl" 

=  P 

(97) 

8£"  +  8£ 

=  7 

(98) 

where  e  is  the  eccentricity  of  the  elliptical  orbit  in  which  the 
guidance  system  is  operating, 

9  is  the  true  anomaly  of  the  vehicle. 


w  r  w  r  w  r 

P  ,  P  ,  and  P„  are  the  coordinates  of  the  error  sources, 

JL  &  v5 

u  =  0,  the  time  rate  of  change  of  true  anomaly,  and 
prime  (*)  signifies  =  i  ^  . 

If,  for  example,  the  accelerometers  of  the  Inertial  Navigation 
System  are  maintained  in  a  local-level  orientation,  then  accelerometer 
bias  corresponds  to  a  constant  input  to  the  equations  of  motion. 
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Figures  (3)  through  (6)  then  may  be  interpreted  as  plots  showing  the 
propagation  of  system  errors  in  nondimensional  altitude  and  cross-track 
due  to  accelerometer  bias,  when 


(99) 


The  sampled-data  solution  discussed  previously  is  of  course  valid  too. 

— 4  — 4 

Hence,  if  a.  «  10  g  ,  a  «  10  g  ,  r  =  100  miles  plus  Earth  radius, 

5  e  H  e 

e  =  0.01,  using  equations  (79)  and  (80)  one  obtains  for  N  =  1  orbit. 


Sx  =  r  S§  **  4  miles 

(100) 

5y  =  r  5*]  ~  -10  miles 

(101) 
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APPENDIX  A.  DERIVATION  OF  THE  TSCHAUNER-HEMPEL  EQUATIONS 


In  this  appendix  the  standard  derivation  of  elliptical  relative 
motion  is  reviewed  for  completeness  and  to  establish  notation.  A 
derivation  of  the  Tschauner-Hempel  equations  j  also  given. 

Consider  the  relative  motion  between  a  reference  object  in  an 

—t 

elliptical  orbit,  described  by  position  vector  R,  and  a  nearby  object 
in  a  slightly  different  orbit,  described  by  position  vector  r.  (See 
Figure  13).  The  relative  position  of  the  second  object  with  respect 
to  the  first  is  designated  by  the  vector  p  so  that 


R  p  =  r 


SaMliic 


/{e/crcncc.  £Uipsc 


Stmriino  Point 


Figure  13.  Coordinate  System  for  Perturbation  Equation 


'l 


II 


r  =  - 


kr 

l?l: 


+  a  (perturbed  body) 


(A3) 


where  k  is  the  gravitational  field  constant,  and  superscript  (  ) 
signifies  differeatietic-n  with  respect  to  time  in  an  inertial  frame. 
Equations  (Al)  and  (A3)  combine  to  form 


a  ♦  g  » -■  + 1 

8  p  if  +  p|3 


(A4) 


By  taking  the  square  root  of  the  dot  product  of  -f  with  itself 
it  is  readily  verified  that 


->  -»  2  -  3/2 

-jj— 3  r„2 2R  •  o  P  vi 

|R  -r  pj  =  [R  (1  + - “■  + 


2  2‘ 
R  R 


(A5) 


2 

If  terms  of  order  are  neglected  as  small  compared  with  terms 


of  order 


/P\ 

l  Ri 


~  —3  -  2?.P  -3/2 

lR  *  PI  =  R  U  + - 5—) 


-3  - 
R  U 


|S  +  ?]-3  =  e'3  (i 


3  / 

—  { - —  )  -r  higher  order  terms  j 


R 


3S  •  p  \ 
E2  7 


(A6) 


(A7) 


With  equation  (A7),  equation  (A4)  becomes 


II  II  ,  ,,  3R’  •  o  ^ 

^  -  ~t  =  -  k(R  +  p)R  (1 - £ — ) 


+  s 


(AS) 
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■yw 


L-TTa-U 


Subtract  (A2)  from  (A8)  to  obtain 


II 

-» 

P  = 


3k  R  .  o 


3  2 

R  R 


k  3k  r  . 

R  "  “3  +  “3  ~ 

R  R  R 


(A9) 


Neglect  the  third  tero  on  the  right-hand  side  as  snail  conpared  with  the 
first  two,  and  the  basic  perturbation  equation  results 


II 

It  -»  3k  -•  — 

p  =  -  r  p  +  — r  (R  •  p  )  R  +  a 

R  R 


(A10) 


L 


If  p  is  the  tine  derivative  of  p  taken  in  the  rotating  reference 

frane  and  w ,  is  the  angular  velocity  vector  of  the  rotating  reference 
wl 

frane  with  respect  to  inertial  space,  then 

I  h 


p  _  p  +  x  P 


(All) 


and 


L 

II  L  I*  L  — >  — >  — >  — >  — >  -4 

?  =  p>+  x  o  *  2Vi  x  p  +  Vi  x(Vi  x 


(A  12) 


In  the  rotating  reference  frane,  if  we  define 


*  x 


A  v 


A  z 


^  R  ’  11  “  R  •  C  R  • 


(A13) 


P  = 


(A14) 


Also  in  the  rotating  frane  we  have 


R  - 


R 

0 


pi  * 


a  = 


i  “t/t 


0 

0 


(A15) 


Combining  equations  (A9)  and  (A12)  ,  resolved  in  the  rotating  frame,  one 
obtains  the  scalar  equations 


(A16) 


CA17) 


(A18) 


where  (*)  signifies  differentiation  with  respect  to  time. 
The  following  identities  can  be  obtained  by  differentiation: 


(A19) 

(A20) 

where  (')  signifies  —  ~  . 

d“  y  dt 

Expressions  identical  in  fora  hold  for  ■'j  and  Combining  these 

equations  with  equations  (A16)  through  (A18) f  and  noting  that 


2'j  e  sin  <r 
1  +  e  cos  r 


(A21) 


yields 


2’lR'C 


z/ri  ’ 


-r  W2R^ 


(A22) 

(A23) 

(A24) 
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Noting  that 


2 

k  fa  R _ 

r2  “  l  +  e  cos 


(A25) 


the  Tschauner-Hempel  equations  are  obtained: 


-  ** 
5 


3  i 


1  +  e  cos 


_  o-  *  —  _ z _  p 

‘  “  2p  1 

v  R 


■  +2^-2lP2 


td2R 


C"  +  U^p3 


(A26) 


(A27) 


(A28) 
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APPENDIX  B.  REVIEW  OF  FLOQUET  THEORY 


In  this  appendix  the  standard  results  of  the  theory  ox  linear  differential 
17  18  19 

equations  *  are  reviewed  for  coapleteness  and  to  establish  the 

notation* 

Theorea:  The  n^-order  linear  inhomogeneous  systea 


x*  (6)  --  F(0>x(0>  +  D(0)u(0)  ;  x(0  )  =  x  (Bl) 

o  o 


has  the  general  solution 


x(6)  =  1(0,  0  )x 
o  o 


X(0,  0  ) 
o 


r 

6 


X^Ct,  0  )D(T)u(T)dT  CB2) 
o 


where  X(0,  0  ),  the  n  x  n  state  transition  natrix,  is  the  solution 
o 

of 


X*  (0.  6  )  —  F(0)X(0,  0  )  ;  X(6  ,  0  )  =  U  (the  unit  catrix)  (B3) 

o  o  o  o 

Proof  1:  Substitute  32  into  Bl. 

Proof  2:  Assume  tfe  solution  x(5)  to  be  cade  up  of  the  complementary 

solution  x  (0)  and  a  particular  solution  x  (0): 

c  p 

x(0)  =  x  (0)  -r  x  (0)  (B4) 

c  p 

The  ntJl-ordered  homogeneous  fora  of  equation  (Bl)  has  n  special  linearly 

independent  solutions  which  can  be  arranged  as  columns  of  an  n  x  n  natrix 

X(0,  6  )  which  satisfies  equation  (B3).  X(0  ,  0  )  was  chosen  as  the 

o  o  o 


*X(0,  0  )  is  known  as  the  catrizant,  fundamental  catrix,  state  transition 
o 

catrix,  or  catrix  of  partials. 
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unit  natrix  so  that  an  arbitrary  cosplenentary  solution  i 


*  -  ». 


x  (S)  =  x<e,  0  )x  (3:.) 

c  o  o 

In  order  to  obtain  the  particlar  solution  assume  the  constants,  >:  ,  of 

o 

the  homogeneous  solution  are  nos  functions  of  9,  and  call  these 
functions  c (0). 


x  (6)  =  XC6,  9  )c (0)  (B6) 

P  o 

This  apparently  arbitrary  assumption  sras  first  cade  by  Lagrange  and 
was  motivated  by  a  desire  to  represent  the  effects  of  planetary 
perturbations  in  the  solar  systea  as  variation  of  the  orbit  elements. 

This  assumption,  it  turns  out,  gives  the  exact  solution  for  the  special 
case  of  linear  equations.  When  (B6)  is  substituted  into  CB1)  we  obtain 

X=c  -t-  Xc*  =  FXc  Du  CB7) 

or 

c*  =  X_1  Du  (B8) 


where  for  an  arbitrary  n  x  n  natrix  A, 


R=0 


(Bll) 


Iesssa:  If  in  the  systea  (Bl),  F(6)  =  F (0  +■  2s),  then  for  any  integer  R, 

X(0  *  2sR,  6  )  =  X(0,  0  *  2s,  S  )  <B12) 

o  o  o  o 

Proof:  X*(0,  9  )  =  F (0)X(0,  9  )  ;  X(0  ,  3  )  =  U  (the  unit  matrix) 

o  o  o  o 

(33).  Since  this  cust  hold  for  all  9, 

Xf(0  +  2s,  6  )  =  F(6  -t  2*}X(0  +  2s,  9  ) 
o  o 

=  F(e)x(e  *  2s,  e  )  CB13) 

o 

since  F(0)  =  F(0  -i-  2s).  The  colunns  of  X(©  2s,  9 >  )  are  n  linearly 

independent  solutions  of  the  hcaogeneous  part  of  (31),  and  therefore, 
each  of  these  columns,  ~  (1  <  R  <  n),  is  given  by  x_  =  X(0,  9  )c_, 
where  for  each  R,  c  is  an  n  x  1  column  natrix  of  constants.  let  C 

A 

be  an  n  x  n  natrix  whose  columns  are  the  cD.  Then 

K 

X(0  -i-  2s,  9  )  =  X(0,  0  )C  (B14) 

o  o 

Since  the  colunns  of  X(0  +  2$,  0q)  are  i-^ependent,  C  ^  exists. 

Equation  (B14)  Bust  hold  for  all  0.  Specifically  it  Bust  hold  for 

0  =  0  : 
c 

X(©  +2 s,  0  )  =  X(0  ,  8  )C  (B15) 

o  o  o  o 

A 

Since  X(0  ,  0  )  =  U.  C  is  known  and 
«>  o 

X(0  -i-  2S,  0Q)  =  X(0,  6o)X(0o  -i-  2s,  0q)  (B16) 


-  52  - 


Equation  (B16)  aust  also  hold  for  all  d.  Specifically  it  eui*1  hold  for 
6  *  6  +  2s: 


X(6  +  Ax,  0  )  —  X(0  +  2k,  0  )X(0  +  2r,  6  ) 

'  /v  '  o  ft  A 


=  x(s,  e  rfie  +  2z,  e  ) 

o  o  o 


By  induction. 


X(0  +  2xR,  0  )  =  X(0,  6  )X*(6  -r  2s,  6  ) 
o  o  o  o 

For  the  balance  of  the  discussion  it  will  be  assumed  that 

F(6)  =  F(6  *  2s). 

Define  a  matrix  »(5,  6  )  by 

o 

Me,  eo)  =  x(e,  t  >€-b(0-0o) 


where  B  is  a  constant  n  x  n  matrix  not  yet  specified. 


x(e,  e  )  =  mb,  e  )€B(e~V 

o  o 


(note  the  similarity  with  equation  (BIO)). 


using  0B12): 


1(6  +  2x,  6  )€B(e+2*_6o)  =  1(6,  0  )cB(6‘eo>  1(6  +  2x,  6  )c 
o  o  o  o 


Bow  define  B  to  be 


B  =  ~  in  X(6  +  2s,  6  ) 

2s  o  o' 


then  from  (B19)  and  (B21): 


1(6  +  2S,  6  )  =  c  *  1(6  +  2s,  6  )c 
o  o  o  o 


1(6  +  2s,  6  )  =  U  (the  unit  matrix) 
o  o 


(BIT) 


CB12) 


(BIB) 


(BIB) 


CR20) 


(B21) 


(B22) 


(B23) 


When  (B23)  is  substituted  into  (B20)  we  obtain  R(0  +  2n,  0  )  =  R(0,  0  ), 

o  o 

so  that  R(0,  0q)  is  a  periodic  mat  -ix.  Now  let 


W  =  eB(6-6o)->w,  _  BW 

(B24) 

Then 

X’  =  FX  implies 

R’W  +  RW’  »  FRW 

(B25) 

R’W  +  RBW  =  FRW 

(B26) 

R*  =  FR  -  RB 

(B27) 

and 

B  =  R*"1  FR  -  R_1R’ 

(B28) 

This  result,  where  F  and  R  are  periodic  and  B  is  constant,  is 
called  rhe  Lyapunov  reduction  of  equation  (Bl). 

Let  A  be  the  Jordan  canonical  fora  of  B;  i.e. 

A  -  "2  CB29) 

then 

w  B<0-0  )  Q~lAQ<©-0  >  .1/  {*-eo>  (B30) 

I  *  C  O  a  c  O  «  Q  C  Q 

and 

X<0,  0  )  -  *<0,  6  )Q'1cA<&“6o)Q  (B31) 

o  o 

If  the  transformation  P(0)  is  introduced  so  that 

z<0)  ,  p'1(0)x(0)  (B32) 


54  - 


where 


z'(0)  =  Az (0)  ;  z (0  )  =  P_1(G  ) x 

o  o  c 


I'.B  »j> 


then 


A(0-e  )  .  A(0-0  )  -1..  , 

z(0)  =  e  o  7/q  )  =  e  o  P  (0  )x 
zv°o  o  o 


(B34] 


Combining  (B32)  and  (B34): 


x(0)  =  P(0)eA(0_eo)  p“1(0  )x 

o  o 


(B/iCv 


From  (035)  it  follows  that 


x(0,  e  )  =  f(0)€a(6"V  ) 

o  o 


(B36) 


If  we  take 


Q  =  P-i (0  )  (B37) 

o 

then  it  follows  from  (B31)  that 

R(0,  6  )  =  P(6)Q 
o 

=  P(0)p"i(6  )  (B3B) 

o 

Equation  (B36)  is  the  form  of  the  state  transition  matrix  used  in  the 
basic  text. 
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f 


a pgagix  c.  derivation  or  a  special  form  of  ras 


SOUfriCn  TO  Z,2*iEAR  8WATJCKS  WHS 

"*■  ■  w»-  I  ■  ’«*  xUo «  ■'»  .«  -  — »■ «•/ 


=£H IODIC  COEFFICIENTS 


The  solution  to 


x‘(^>  -  «•(£)- 


C5) 


r  ^  ;  .  /.  % 


s 

r. 


has  >>«?«*  sbo^vii  \o  b*  (see  ^poc.-'.^ix  S> 


\(3)  « 


'J  D 

c  o 


,1  \»*  , 


..9 

|  *  ^(t,  9uJ>D(t)u(t)cIt 


r 


O 


NT  - 


*o 


V- 

* 


— *— 
<*'r 


A 


— V 


— ■ 


r 


4^+MT 


♦ . -+ - ■■■» 

♦ 

'■» 


* 


Lcm«:  If  F(3)  =  Ffe  +  T)  then 

Xfe,  6  )  =  X(£>  +  a,  e)X?(S  +  T,  s  ) 

o  o  o  o  o 

Proof:  In  Appendix  B  it  vas  established  (for  T  =  2s,  no 

that 


xfe  +  st,  e  )  =  x<9,  e  )***&  +  t,  e  ) 


Let 


S  =  S  +  XT 


(Cl> 

(C2) 

(C3) 


(C4) 

restriction) 

(C5) 

(06) 
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iher: 


X(S 


-  X(S 


A  c- 


r 


v*  ; 


Fross  <C3) 


x&.  ;  >  -  X(t  ^  j,  5  »:*  Cv  ♦  T,  £  > 

n  o  o  o  o 


(Cf  } 


Sot;  tj  *:'*iag  rc!«tiioj.  v'lb  t  i  ^cj iio-i  ^C2/  yields. 


r  /-  t*w  ^ 


:*.  5  )**(€- .  -*■  T,  g  »X, 


— >\- 


/  “Ji'ri :'•**  JC9) 


Leaaa:  lx  flfe)  =  D(0  +  7)  and  u(0)  **  ud6  +  T>  the  solution  (C9)  can 


be  written 


x(S)  =  X{9  +  c,  0  )C*jr  +  Xfe  ♦  o,  0  > 

o  o  o  o  o 


N  x 

k=l 


+  X(s  +  o,  0  ) 

o  c 


/•v° 

J 


e  >D(T)u(T)dT  (CIO) 

o 


whers 


C  =  X(0  +  T,  0  ) 

o  o 


(Cll) 


*1  “ 


9  ■‘■T 

J'-' 


X  ^ (t,  0  }D(T)u(x)dT 
o 


(C12) 
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**  **♦*  1 


Proo: 


S  '(• 


e_ 


5  }  D(t)u  (  "Odl 

r 


9  +T 

o 


J 


X  iCx,  S  /.  (t)u(t)c1t 
o 


5 


4.51 


-1 


f 


X  (t,  C  (-r)dT 

v» 


1 1 . 


§  *:<T 
o 


/  +<K-1}7 
o 


X  fx  )D<T)u(T>dT 
o 


e 


J 


X_1(T, 


e  +JTT 
o 


e  )D(t)u(t)c!t  . 
o 


By  sjjiple  changes  of  variable  in  each  integral  obtain 


X_1 (t,  V  )D(T)uC')d- 
o 


J 


$  +t 

o 

X 


-1 


e 


o 


(T, 


0  )D(t)u(t)<It  + 

c 


(see  page  59> 


58 


u 


Proof : 


|  X  “(t.  5r>C(x)«< '.)dt 


Q  +T 
t  ° 

J  X  ^T, 
$ 

o 


6  /.  (t)u(T>dT 


8^27 

j"  X  1(t,  C^iZwiV*. T/dt  + 

9. 


^  «  »  • 


y  o 

j  X“  (ts  f>  >D(t)u(T>dT 

V  ° 

fr  +(S-1}7 
O 


8 

+  J"  X  1(t,  8o)D(x)u(T)dT  . 

e  +nt 


By  staple  changes  of  variable  In  each  integral  obtain 


X-1(t,  6  )D(x)u(T)d- 
o 


J 


$  +T 
o 


-1 


e 

o 


(T, 


8  )D(x)u(T)dT  ♦ 

O 


(see  page  59> 


58  - 


e  +t 

JVl 


(t  +  T,  $  ) D( i )u ( " )  —  t  ^ 
o 


6  +T 
o 


f  -1 

..  +  j  X  (x  +  (X  -  1)T,  6o)D(t)a(T)d: 


e 


e-a n 


j  X_1(t  +  XT, 


9  )D(x)u(x)dx 
o 


9 


Use  X(t  +  XT,  9  )  =  X(£},  3  )X*  (g  +  T,  Q  )  shown  above  to  obtain 
o  o  o  o 

X_1 (x  +  XT,  9  )  =  X~*(S  +  T,  9  )X_i (x,  9  ) 

o  o  o  o 


o 

8, 

9_ 


9  +T 
o 


>D(x)u(T)d 


-  J 


9  )D(x)u(x)dx 
c 


9. 


9  vr 


+  X_1te  +  T,  9  ) 
o  o 


/  °  x-c. 


9  }D(x>*  (x)dx  + 
o 


9  +T 

♦ - +  x"<Jf“1)fe  +  T,  $  )  f  X"1^, 

o  °  J 


9  )D(r)u(x)dx 
o 


9  +c 
o 


-k  -i 

X  (9q  +  T,  9o)  j  X  (x,  9o)D(x>u(x)dx 
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Introduce  this  last  relationship  into  equation  (C9)  and  simplify 
to  obtain  expressions  (CIO)  through  (C12).  Thus  are  proven  expressions 
(36)  and  (37)  of  the  basic  text. 

Leana:  If  D(0)  =  D(g  +  T)  and  u(9)  =  u(9  +  1CT)  where  M  is  an 

integer  the  solution  (equation  (C9))  can  be  written 


xfe)  =  x(s  +  o. 


><*, 


9  x 
o  o 


+  xfe  +  o, 


e  ) 

o 


r-1 

(£ 

k=0 


Q-rXI 

+  Xfe  *  c,  9  )CN"rM  (  X_i(T,  9  )D(x)u(x)dx  (C13) 

o  o  J  o 

9 


where  C  =  X(S  +  T,  9  )  (C14) 

o  o 

e  4«r 

V  =  I  X_1(T,  8  )D(t) u(x)dx  (C15) 

2  J  o 

o 

r  is  .  6  in  eger  such  that  rM  <  N  <  (r  +  1)11  (C16) 

6 

(x,  9  )D(x>i(x)dx 
c 

o 


Proof: 


f-1 


e  +*r 


j 


-i 


(x,  9  )D(x)u(x)dx  * 
o 


S^+rKT  9 

...  +  f  x”X(t,  e  )D(x)u(x)dx  +  f  X~" 1  ( T .  0  )D(x)u(x)dx  . 

J  °  Jo 

5  +(r-l  )MT  g  +riCr 

c*  o 


-  60  - 


i 


where  r  is  described  by  equation  (C16)  and  X  lv  'till 
largest  number  of  integer  values  of  T  in  6  .  Again  u.-t  .-lnple 
variable  changes  in  the  integrals  to  obtain 


0 

J*  X  *(',  0Q)D(T)u(T)dT  = 

e 

o 


6  +MT 


a 

~  o 


d  )P(T)u(:)d-  + 
o 


6  +KT 

f*'1 


(t  +  (r  -  1)JC,  0o)D(T)u(x)dT  + 


6 


0-rlff 


/  X-1(t  +  riff. 


9 


6  )D(t)u(t)c!t 
o 


Again  X  *(t  +  W,  6  )  =  X  +  T ,  0  )X  ^(t,  6  ),  so  that 

o  o  o  o 


6 


e  +xr 

o 


J  X-1(x,  0o)D(x)u(')dx  =  j  X_1(t, 


6  )D(-r)u(t)dx 
o 


6 


*  x‘M( e  t,  0  ) 

o  o 


0  +MT 

•  c 


J  X'^T, 

9 


6  )D(x)u(T)dT  + 
o 


+  ....  x"(r"i)MC 0  M,8) 

o  o 


(see  page  62) 


0  +MT 
o 


/« 

x 


9  )D(T)u(x)dT  4* 
o 


e 
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‘■■ST*  *’,“!«* 


-rSJ 

+  X  (9  +  T,  9  ) 

o  o 


0-rJfT 


J  X_1(T, 


0o)D(-)u(T)dT  . 


Substitution  into  equation  (C9)  and  siaplif ication  yields  relations 
(C13)  through  (C15).  Thus  equations  (39)  and  (40)  of  the  basic 


text  are  proven. 

t:  If  Dfe)  =  Dfe  +  T)  and  u(0)  =  u (9  +  P)  where  P  t  MT  for 


M  =  0,  lf  2,  . .  then  define  K  such  that  KP  =  JfT,  where  K 

and  M  are  both  integers.  Then  x(0)  nay  be  approximated  by: 


x(9)  =  X(S  +  c,  0  )C*Nx  +  X(9  +  a,  0  )  (  V* 

o  co  o  o  \  J 


0-rKP 


X(g  ♦  a,  0  )CN_rM 
o  o 


J  x^Ct. 


0  )D(')u(t)dT 
o 


(C17) 


where  C  =  X(0  +  T,  6  ) 

o  o 


KP  =  KT,  rKP  <  NT  <  (r  -  1)KP 


(C19) 


0  +KP 
/•o 


■  /  ^ 


0  )D(-r)u(t)d' 
o 


(C20) 
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Proof:  The  proof  of  relationships  (C17)  through  (cj  » >  :  .. 
the  proof  for  relationships  (C13)  through  (C16)  directlj. 


\ 
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APPENDIX  D.  THE  I -MATRIX 


i 


The  samoled-data  solution  defined  by  equation  (73),  that  is 


where  p  is  the  seuilatisrectum  of  the  reference  (nominal)  ellipse 
k  is  the  gravitational  field  constant,  and 
e  is  the  eccentricity  of  the  reference  ellipse. 


For  disturbances  of  the  form 


u(0)  = 


j\9 


jK2e 


where  K  and  K  are  integers, 

JL  4 

c1  and  c  are  complex  constants, 

X  4d 

the  I-matrix  becomes 


_  (2  +  e) 

(1  +  e)f1  q  +  e)  f2 


/-.  \  2^3/2 

(1  -  e)(l  •  c  ) 


(3ef3  "  *4^ 


2  2 
e(l  -  e)  1  -  e 


{3ef 3  -  [1  -  (1  -  e)2 \l  -  ?  ]  f J 


<2  +  e)  fa  +  -jf-f-fj  f2 


where  p,  k,  and  e  are  as  defined  in  the  previous  case,  and 


16e(l 


^  {j32c1(l-.2)2K1Z*1  -  ,3c2Z*2  4C«o4«2)Z* 

-e  )  V 


+  «^-K2)Z l  -  2e(7K2-2)Z®  -  4<K^+9K2-12)Z*  -  2e(l^-2<»Z* 
-  4(I^«9K2-12)Z2  +  2e(7K2+2)Z2  +  4(^+1^^  +  e<l^-X2)]j 
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/  K  K 

=  ^2  { jClKlZl  “  °2Z1  } 

=  -  "'i  a  j2*2-!2!1 2<Kizi‘4ZrKi) +  J°2k2zi2  2t«VJ 

4e(l-e  )  l 

+  4e(K2-l)Z^  +  4<l-2e2)Z2  -  4e(K2+l)Z]L  -  e2K2]J 


/  K  K 

|c1Z11(l+K1ln|Z1{+jKlIr)  +  c1K1S1  +  c^CjlnlzJ-*) 

+  ^VVJ  +  XT  (C1  +  -  C2S2<V3} 

'  2V5/2  f'Vl1  4  (t(8e+4'3)  ^2H«^-3K1+2)Z® 


2.  2x5/2  ‘ 

48e  <l-e  )  V. 


-  (2K^-8)Z*  +  (K^+3K1+2)Z^]  +  [2  -  2(l-e2)3/2  HO^-lC^zJ 

-  (K2+9K1-12)Z^  -  (^-9Kr12)Z3  +  (K^+K^ZJ  + 

+  [e  -  (e~2e3)  ’Jl-?]  [(K^+K^Z®  -  <14X^4  )Z®  -  <2K^-40)Z* 
+  (14K1+4)Z2  +  (K^-K1>  3  +  [24e(l-e2)1  [(K^JZ® 


(DIO) 


-  ©K^Z*  -  K.Z^  -  j4c2(l-e2)K2Z12  ^[(4+2e2)  1 1 

-  4]  [e(K2-l)Z3  -  c(K2+1)Z1]  +  [(l+2e2)  -1]  [e^Z* 

+  <4-8e2)Z^  -  e2^])  j 
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r-1 


-ETPiofe)  {zi«-vr'k  -  \k<-vr"k  -z2a-vr'k 


k=0 


_k, 
+  Z2( 


(K1) 


uv 

(2k)  I _ /  jeV 

“  22k(k:)2(2k+i)  ' 2  ' 


2k+l 


k=0 


/Li_ 

^j(2k+l): 


K  -2k+l 

x  d2k+l  (  (Z2-l)2k+1(ZiZ2-2Z+Z1)2k+1Z  1 


(2k+l):  dZ2k+l  \ 


2k+2,_  _  v2k+2 

e  (Z-Z2) 


'  Zk-k  . 

1  d _ _  / 

(2k-K  ) l  2k-r  l 
**  dZ 


I 


z=z. 


2k"Kl  /  (Z2-l)2k+1(Z  Z2  -2Z+Z  )2k+1 


(eZ2  +2Z4e)2k+2  /j^  J 

KjS2k 


z  =  -  —  (1  -  Ji  -  «2)  , 

1  e 


z 2  =  “  "T"  (1  +^X  “  e2)  • 


The  infinite  series  Sg  in  the  f3  ten  above  arises  in  the 


tlon  of 


r  2s  __-l 


„  1  /  sin  x  jKx. 

VK)=— j0  <i  7.  So.  t>  £  dT 


x.  = 


e  sin 


(1  -  Zj  cos  t) 


(Dll) 


(D12) 

(D13) 

(D14) 

evalua- 

(D15) 


where 


1  +  e  cos  t 


(D16) 


An  approximation  to  this  integral  can  be  made  by  observing  that 


X  -  e  sinr  [l 


2  2  3 

(Zj+e)  cost  +  eCZ^e)  cos  x  -  e  (Z^+o)cos  t  +  . ...]  (D17) 


and 


-1  13 

»ln  X  =  X  +  g  X  + 


+ 


(2k): 

2k  2 

2  <kl)  (2k+l) 


2k+l 

X  +  .  .  .  4 


(D18) 


Since  Z^  and  e  are  of  the  sane  order  of  magnitude,  x  ■*>’  be 

approximated  to  whatever  accuracy  desired  by  cutting  off  the  series  and 

3  5 

discarding  similar  powers  of  e  in  x  »  X  *  etc.  Term  by  term  integration 
can  then  be  accomplished  on  the  unit  circle. 
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